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Abstract. A model of strongly inhomogeneous medium with simultaneous perturbation of 
t"-- ■ rigidity and mass density is studied. The medium has strongly contrasting physical character- 

C> istics in two parts with the ratio of rigidities being proportional to a small parameter s. Addi- 

tionally, the ratio of mass densities is of order e _1 . We investigate the asymptotic behaviour 
of spectrum and eigensubspaccs as e — > 0. Complete asymptotic expansions of eigenvalues and 
cigenfunctions are constructed and justified. 

We show that the limit operator is nonself-adjoint in general and possesses two-dimensional 
Jordan cells in spite of the singular perturbed problem is associated with a sclf-adjoint operator 
' in appropriated Hilbert space C e . This may happen if the metric in which the problem is 

■ self-adjoint depends on small parameter e in a singular way. In particular, it leads to a loss of 

p I ' completeness for the eigenfunction collection. We describe how root spaces of the limit operator 

C/_ ■ approximate eigenspaces of the perturbed operator. 

•i— > . 

c3 . Introduction 
g- 

l_j | We consider a model of strongly inhomogeneous medium consisting of two nearly homogeneous 
components. Assuming a strong contrast of the corresponding stiffness coefficients k\ <C k 2 , we 
get that their ratio ki/k 2 has a small order, which we denote by e. In general, the mass densities 
ri and r 2 in two parts could be quite different as well or could be the same. We model this 
assuming that the density ratio r\jr 2 is proportional to E~ m . We investigate how the resonance 
vibrations of the medium change if the parameter e tends to 0. In the one-dimensional case we 

£sj . consider the spectral problem 

r " ' d ( du \ 

— lk £ (x)-j^-j + \ £ r £ (x)u £ = in (a, b), aiu' e {a) + a u £ (a) = 0, Piu' e (b) + (3 u £ (b) = 0, 

• rH ! where (a, b) is an interval in R containing the origin and 
^ ' ( ( 

H! , / s k(x) for x G (a, 0) . . e~ m r(x) for x G (a, 0) 

- ' k ^ = \ \ \ f c n h\ r ^ = \ n f c n h\ ^ 

\eh{x) torxG(U,o), p[x) tor xG (U,o). 

Here k, r and x, p are smooth positive functions in intervals [a, 0] and [0, 6] respectively. 
At point x = of discontinuity of the coefficients we assume that transmission conditions 
w £ (-0) = u e (+0), (ku' E )(-0) = £(xu' e )(+0) hold. 

Of course, the limit properties of spectrum depend on the power m characterizing the density 
ratio. Intuitively, we expect that for large values of m the mass density perturbation has to be 
dominating whereas for small m the rigidity perturbation has to be leading. Then it has to be 
at least one critical point m separating the cases. It appears to be truth exactly for m = 1, 
when the mass density perturbation is strictly inverse to the stiffness one. 
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This paper is devoted to the critical case m — 1. We consider the Dirichlet problem 



(k(x)u' e )' + e~ 1 X £ r(x)u £ = 0, x E (a,0), (2) 

e(x{x)u' e )' + X £ p{x)u £ = 0, xe(0,b), (3) 

Ue(-O) = u E (+0), (ku' £ )(-0) = e (x<)(+0), (4) 

u £ {a) = 0, w £ (6) = (5) 
and investigate the asymptotic behavior of eigenvalues A e and eigenfunctions u e as e — > 0. 



After a proper change of spectral parameter problem ©-© can be represented as a problem 
with perturbation of the transmission conditions only (cf. the example with constant coefficients 
below). At first blush, the problem looks very simple. But the point is that the problem shows a 
complicated picture of the eigenspace bifurcation. In Section [2] we prove that the limit behavior 
of the spectrum is described in terms of a nonself-adjoint operator that has in general multiple 
eigenvalues and two-dimensional root spaces. At the same time, (J2])-© is associated with a 
self-adjoint operator in the weighted space C e with the following scalar product and norm 

O,V0e = + (P0> ^)ia(0,6), \\4>\\e = VWWe ■ (6) 

It is obvious that for each fixed e > the spectrum of ([2])-([5]) is real, discrete and simple, 
< Af < Af < • ■ ■ < \j < ■ ■ ■ — > oo as j ' ^ oo and the corresponding real- valued eigenfunctions 
{u E j}°^ l form an orthogonal basis in C e . How may it happen? The metric in C e for which 
the perturbed problem is self-adjoint, depends on small parameter e in a singular way. In 
Sections El S]we construct and justify the complete asymptotic expansions of eigenvalues and 
eigenfunctions. Therefore there exist pairs of closely adjacent eigenvalues A^ and A^ +1 being 
the bifurcation of double limit eigenvalues. Although the corresponding eigenfunctions u e j and 
Mej+1 remain orthogonal in C £ for all e > 0, they make an infinitely small angle between them 
in L2(a,b) with the standard metric and stick together at the limit. In particular, it leads to 
the loss of completeness in £12(0,, b) for the limit eigenfunction collection. Nevertheless both u £ j 
and u e j + i converge to the same limit, a plane ir(e) being the linear span of these eigenfunctions 
has regular asymptotic behaviour as e — > 0. In fact, a root space 7r corresponding the double 
eigenvalue is the limit position of plane ir(e) as e — * 0, as is shown in Theorem We actually 
prove that the completeness property of the perturbed eigenfunction collection passes into the 
completeness of eigenfunctions and adjoined functions of the limit nonself-adjoint operator. 

This work was motivated by [H Ch.8], where the similar problem for the Laplace operator 
has been considered. The authors have handled the limit operator as the direct sum of two 
self-adjoint operators that nevertheless does not entirely explain the bifurcation picture in 
perturbation theory of operators. The aim of this paper is to present more rigorous and 
detailed study of the case in operator framework. 

Finally, let us remark that the vibrating systems with singularly perturbed stiffness and 
mass density have been considered in many papers. In the case of purely stiff models (with 
homogeneous mass density), the asymptotic behavior of spectra have been studied in [B] - [12J. 
Referring to problems with purely density perturbation often involving domain perturbations, 
we mention [13]- [18] with the latter including a broad literature overview in the area. Spectral 
properties of vibrating systems with mass entirely neglected in a subdomain were also studied 
in [19], [20]. To the best of our knowledge, the first asymptotic results for the problems with 
simultaneous perturbations of mass density and stiffness appear in [21], [22]. 
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1. Preliminaries 

We demonstrate an example where eigenvalue bifurcation is calculated explicitly. If all coef- 
ficients in (J2]), ([3]) are constant we get the eigenvalue problem 

y'; + cu 2 e y £ = 0, i6 (a,0)U(0,6), (7) 
y £ (a)=0, y E {b) = 0, y e {-0) = y £ {+0), y' £ (-0) = ey' £ {+0), (8) 

where u £ = £ _1 A e . Then each non-zero solution can be represented by 



Ve 



A £ sin uj £ (x — a) foi x E (a, 0), 
B £ $mu) £ (x — b) for x E (0, 6), 



with oj £ > and A £ ,B £ E R. By virtue of (jHJ) we have 

A e sin oj e a — B £ sin u e b = and A e coso; e a — £.B e cosu E b = 0. 

Looking for a non-zero solution of the algebraic system, we get the characteristic equation 

coso; e a sin u £ b = esina; e a cosu £ b. (9) 

The latter easily gives existence of the limit uj £ — > uj as e — > such that 

cosa>a sincufr = 0. (10) 

Moreover, the root uj has to be positive. Obviously, if we suppose, contrary to our claim, that 
uj £ goes to as e — > 0, then ([9]) can be written in the equivalent form 

cosa; e a smuj £ b 



cos u £ b sin uj £ a 



for sufficiently small e. A passage to the limit as e — > and uj £ — > leads to a contradiction, 
because the left-hand side converges towards the negative number 6/a. 

If a and 6 are incommensurable number, then all roots of (jlOp are simple. In fact, multiple 
roots exist iff 2n|a| = (21 — l)b for certain natural / and n. Let us consider the case a = — 1 and 
6 = 2. Then the lowest positive root uj = n/2 of (jlOp has multiplicity 2. On the other hand, 

equation (j5J) admits the factorization ^cosu; e — ^2+2^) ^ COSCo, e + x/irS') smu; e = 0- Hence 
the lowest eigenvalues u e> i = | — arcsin -v/^r^, ^ e ,2 = § + arcsin ^2+2! are c l° se ly adjacent 
and converge to the same limit 7r/2. The corresponding eigenfunctions and y £i 2 are defined 
up to a constant factor as 



(x ) = J (-i)V 2e /( 1 + £ ) sin ^-( a; + 1 ) for xe (-1,0), 
yejW 1 anw ej -(a;-2) for x 6 (0,2). 1 J 

We see at once that the angle in L 2 (— 1, 2) between the eigenfunctions y £> \ and y Ej 2 is infinitely 
small as e — > 0, because both eigenfunctions converge towards the same function 

f for XE (-1,0), 

~ isin-(x-2) for xE (0,2). 
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The point of the example is that the collection of eigenfunctions {u e j}'jL 1 loses the com- 
pleteness property at the limit on account of the double eigenvalues. We now turn to per- 
turbed problem (E])-([5]) in the general case. To shorten formulas below, we introduce notation 
I a = (a,0), h = (0,6) and 

K{x) = [ k{ ?\ iOIXeIa R(x) = { r{ *\ iOTXeL 
I k[x) for x G I b , I p[x) for x G h- 

Proposition 1. For each number j G N eigenvalue A^ of is a continuous function of 

e G (0, 1) and ce < A^ < CjE with constants c, Cj being independent of e. 

Proof. The continuity of eigenvalues with respect to the small parameter follows immediately 
from the mini-max principle 

f° kv 12 dx + e L kv' 2 dx , 
X] = min max — 5 ^ , (12) 

Ej vGEj £ -l f rv 2 dx -{- f OV 2 dx 
v^O Ja J ° r 

where the minimum is taken over all the subspaces Ej C H^(a, b) with dim.&, = j. We consider 
the eigenfunctions v\, . . . , Vj corresponding to the lowest eigenvalues pi, . . . , pj of the problem 

{k(x)v')' + pp(x)v = 0, x G I b , v{0) = v(b) = 0. (13) 

Extending each Vk by zero to (a, 0) we get that the span M. of V\, . . . , Vj is an j-dimensional 
subspace of Hl(a,b). Then 

f° kv' 2 dx + e f„ 6 kv 12 dx e f„ 6 xv' 2 dx , 

X" < max -is — ^o = ma x = ep h 14 

' «eA4 £ -i J\ v 2 dx + j h pv 2 dx ^ M j Q pv 2 dx 

which establishes the upper estimate. Next, by the same mini-max principle 

f° kv' 2 dx + e C kv' 2 dx k* f ° v' 2 dx + e f„ v' 2 dx eK uo 2 , 
A^ > \\ = min — 5 ^_ > _ m i n — _ £^ > ^ 

ffo( a > 6 ) J a ry 2 + J pv 2 dx r * H o( a > b ) e~ l f v 2 dx + f Q v 2 dx r * 

where fc* = mm x ^ a ^ K (x) , r* = max x6 ( a b ) i?(x) and u;^ is the first eigenvalue of problem 
(J7|)-(jSD with constant coefficients. It remains to note that u Ej % — > tt/2. □ 

2. Convergence Results and Properties of Limit Problem 
Let us consider the eigenvalue problem 

(K(x)u')' + pR(x)u = 0, x G I a U I 6 , 

u(a)=0, «(6) = 0, u(-0)=«(+0), u'(-0) = 0, 1 J 

that will be referred to as the limit spectral problem. The spectrum of ( TToT) is discrete and 
real (see Th. [TJ below). We introduce the space Ti = {/ G Hl(a,b): f a G H 2 (a, 0) and /& G 
if 2 (0,6)}, where / a and /& are the restrictions of / to intervals I a and resp. Problem (fTol) 
admits the variational formulation: to /ind /i G C and a nontrivial u 6?i sitc/i £/ia£ 

/•& /•& 

/ Ku'<p'dx + k(0)u'(+0)<P(0) = p Ru<pdx (16) 
/or a// G C^°(a, 6). We first prove a conditional results. 
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Proposition 2. Given eigenvalue X s and the corresponding eigenfunction u £ of (0|)-([3J), if 

e~ l \ £ — > fi* and u £ — > in H 2 weakly on each intervals I a , I b and u* is different from zero, 
then fi* is an eigenvalue of fT5\) with the eigenfunction u* . 

Proof. We make a change of spectral parameter X £ = Efi e in ([2])- whereat we can reduce 
equation ([3]) by the first order of e. Then each pair (//, u e ) satisfies the integral identity 

rb rb 

/ Ku' e <j)'dx + (1 -£)x(0)<(+0)(/.(0) = // £ / Ru £ (j)dx (17) 

J a J a 

for all <j) G C^°(a, b). The weak convergence of u £ in H 2 (0, b) gives the convergence u £ — > in 
C 1 (0, 6), in particular, u' e (+0) — > w^(+0) as well as 0) — > 0. Moreover, the limit function 
belongs to H., since each u £ is a continuous function at x — 0. A passage to the limit in 
(]17p implies that pair (p*, it*) satisfies identity (T3~6]) . Recall that it* is different from zero, which 
completes the proof. □ 

Before improving the convergent results, we first compute the spectrum of the limit problem. 
Let us introduce space £ = L 2 (r, I a ) ©L 2 (p, h), where L 2 (g, I) is a weighted L 2 -space with the 

1 /2 

norm ||u|| = (Jj g\v\ 2 ) . We consider two operators 

A i = ~ l r£ k £ in L 2(r, Q, V{A X ) = {ue H 2 {I a ) : u(a) = 0, u'(0) = 0}, 

A * = - l A*± in L 2(p, h), V(A 2 ) = {ue H 2 (I b ) : u(b) = 0}. 



— IA. 


^ dx 


p dx 


y (USD 


we f 







(t 


A 2 



-4= (7 m£, D(A) = {Kuj) £ W ® D(4): m(0) =«2(0)}. 

The operator A is nonself-adjoint. Actually, it is easy to check that 

A* = jQ , V(A*) = {( Vl ,v 2 ) g V{A 1 )®V{A 2 ): (A^)(0) = (*t£)(0)}, 

where A\ is the extension of operator A\ to X'(Ai) = [u G if 2 (a, 0) : it(a) = 0} and A 2 is the 
restriction of A 2 to T>(A 2 ) = {m6 T>(A 2 ): it(0) = 0}. Let cr(A) and g(A) denote the spectrum 
and the resolvent set of an operator A respectively. Let IZ^A) denote the resolvent (A — fil)~ l 
of an operator A, where X is the identity operator in C 

Definition. Let u be an eigenvector of A with eigenvalue fi. A solution to {A. — fi1)u* = u 
is called an adjoined vector of A (corresponding to the eigenvalue fi). 

Theorem 1. (i) a (A) = (t{A x ) U a(A 2 ). 

(ii) If fi belongs to cr(A) \ (c(Al) r\a(A 2 )), then fi is a simple eigenvalue. If fi G cr(Ai) C\a{A 2 ), 
then fi has multiplicity 2 and the corresponding root space is generated by an eigenvector and 
an adjoined vector of A. 

(Hi) The set of eigenvectors and adjoined vectors of A forms a complete system in C 

Proof, (i) Let us consider the equation (A — fiX)u = f for fixed / G C. In the coordinate 
representation we have A\ u\ — fiU\ = fi, A 2 u 2 — fiu 2 = f 2 . If fi G" cr(Ai), then u\ = IZ^A^fi. 
In order to find u 2 we introduce the bounded intertwining operator T M : H 2 (I a ) —* H 2 (Ib) that 
solves the problem (xip')' + fipip = in I b , ?/;(0) = g(0), ip{b) = for each g G H 2 (I a ). Note 
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that T M is a well-defined operator for all /i G q(A 2 ). Then u 2 = T^TZ^Ai) fi + TZ^(A 2 )f 2 and 
the resolvent of ^4 can be written in the form 

nM) = L „ ?, J • d8) 



T^(Ax) 7^(A 2J 

From the explicit representation of IZ^A) it follows that sets <j(A) and cr(Ai) U<t(t4 2 ) coincide. 

(ii) We suppose that p G o~(Ai) \ cr(A 2 ). Then there exists an eigenvector = (iti, T^iti), 
where u\ is an eigenvector of A\ and, that is the same, one is an eigenfunction of problem 
(k<py + jircf) = in I a , 0(a) = </>'(0) = 0- Note that /i is a simple eigenvalue of the problem. 
Indeed, (A — pI)Ufj, = follows from the evident equality (A 2 — /iX)T M = for all p G g(A 2 ). 

Suppose now that p G c(A 2 ) \ cr(Ai). Then operator A has the eigenvector = (0,u 2 ), 
where w 2 is an eigenvector of A 2 . In other words, u 2 is an eigenfunction of the Dirichlet problem 
( TT3l) . Note that each point of u{A 2 ) is a simple eigenvalue. Furthermore, the first component 
ui must be zero, since fi G" cr(Ai). 

Finally we shall show that each point of intersection a(Ai)r\a(A 2 ) is an eigenvalue of algebraic 
multiplicity 2. Obviously, vector = (0,u 2 ), which appears above, is an eigenvector of A in 
this case too. Next we consider the system 

Aii>i - jj vi = 0, A 2 v 2 — /i v 2 — u 2 (19) 

determining adjoined vectors. If vi = 0, then v 2 must be a solution of the boundary value 
problem {>«$')' + fip<f) = —pu 2 in 1^, 0(0) = <p(b) = 0, which is unsolvable. Actually, since 
/i G a(A 2 ), by the Fredholm alternative the problem admits a solution iff J Q p|w 2 | 2 dx = 0. This 

contradicts the fact that u 2 is an eigenvector of A 2 . Consequently we have to assume that v\ 
is an eigenvector of Ai and examine the problem (>cv 2 )' + ppv 2 = —pu 2 in lb, v 2 (0) = fi(0), 
v 2 (b) = 0. Here the Fredholm alternative gives the solvability condition 



(0)4(0)^(0) = - / pu\dx. (20) 



o 



We satisfy one by normalization of v±, because u 2 (0) is different from zero. This condition 
assures the existence of v 2 and a solution V* = (vi,v 2 ) of system f[T9"j) . Vector V* is the 
adjoined vector of A. Pair {V^, V*} forms a basis in the root space that corresponds to p. 
The last statement of the theorem follows from the Keldysh theorem [3] . □ 

We investigate the limit behaviour of eigenfunctions u £ ^ n normalized by conditions 

R{x) u 2 s>j (x) dx = 1, u' £>j (b) > 0. (21) 

Let us enumerate the eigenvalues of operator A in increasing order and repeat each eigenvalue 
according to its multiplicity: p\ < p 2 < ■ ■ ■ < pj < ■ ■ ■ . The next statement improves the 
conditional results of Proposition [2J 

Theorem 2. There exists a one-to-one correspondence between the set of eigenvalues {\ £ j}°^ 1 
of perturbed problem (J2j)-(j5]) and the spectrum of operator A. Namely, £~ x Aj — > pj as e — > 0, 
for each j G N. Furthermore, a sequence of the corresponding eigenfunctions u £ j converges in 
H 1 (a,b) towards the eigenfunction u with eigenvalue pj. 
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'22) 



Proof. For the perturbed problem (j2j)-^5]) we assign the matrix operator in C 

Ae= { A X)' P (^) = {(^' n 2)G^(ii)©^(A 2 ): 

Ul (0) = u 2 (0), (K)(0) = s(W 2 )(0)}. 

Clearly, if /x e belongs to cr(.A e ), then e/i £ is an eigenvalue of ©-(jSJ). Let us solve the equation 
(A £ — jiX)u — f for / = (/i,/2) G £ and // G £>(*4. £ ). Similarly to the previous theorem 
we obtain u x = K^A^fx + eS^, ^2 = + TZ^(A 2 )f 2 , where S^. H 2 (I b ) -> # 2 (/ a ) is a 
bounded intertwining operator that solves the problem (kip')' + /zr?/> = in 7 a , ?/>(a) = and 
(ki[)')(0) = (xg')(0) for each g e H 2 (I b ). This yields that 

/ -eSA f Ul \ (n,{Ai)fi 

-T, I ) \u 2 ) [n,(A 2 )f 2/ 

where the matrix operator in the left-hand side is invertible as a small perturbation of the 
invertible one. Letting e — > we can assert that 

" [ e) \-T, I J \ TIM,)) \T, I J { K»{A 2 ) 

Hence, 1Z^{A £ ) —> 1Z^(A) in the uniform operator topology as e — > 0, which establishes a 
number-by-number convergence of the corresponding eigenvalues [31 Th. 3.1]. 

Next we prove existence of the limit for the eigenfunctions under normalization condition 
(I2TT). We conclude from (ED that f a K(x)u'*(x) dx + (1 - e)x(0)<(+0)n e (+0) = /i e . For each 
v there exists a twice differentiate solution ip(x,u) of equation (xv')' + ispv = in I b that 
satisfies conditions v(b) = 0, ?/(&) = 1. Moreover, ip(x, u) is an analytic function with respect to 
the second argument for each fixed x (2j Th.1.5]. In particular, tf)(x, pf) — »■ ip(x, fi) in C 2 (0, 6) as 
fi e — > yU. Then there exits constant /? e such that u £ (x) = j3 £ ip(x, fi £ ). Moreover, (3 £ is bounded as 
e — > 0, which is due to condition (I2ip . Therefore the values u £ (+0) and m^.(+0) are bounded with 
respect to e. Consequently we have J b K(x)u' 2 (x) dx < /j, e + (1 — e)x(0)\u' e (+0)u e (+0)\ < M. 
Then finally the sequence {u £ } £>0 is precompact in the weak topology of if 1 (a, b). Let us 
consider a subsequence it E / such that u £ > — > u in H l (a, b) weakly. We get u e '(x) = (3 e /if)(x, // ) — > 
/3ip(x,n) = u(x) in C 2 (0,fc) for certain /?. Note that /3 > 0, which is due to fl2~lj) . Moreover, 
Mg/(+0) — > «'(+0) as e' — > 0. A passage to the limit in (|T7|) implies that partial weak limit m 
satisfies the identity 

rb rb 

/ K(x)u'<t/dx + x(0)u , (+0)<l>(0)=[i R{x)u(j)dx 

J a J a 

for all G Co°(a, 6). Moreover, it is different from zero, since J a 6 -R|w| 2 c/x = 1. Consequently 
each weakly convergent subsequence of {u £ } £>0 tends to u, where u is an eigenfunction of ffT5l) 
that corresponds to the eigenvalue \i and satisfies conditions ||«||i2(_ R ^ a ^\ = 1 and w'(6) > 0. 
Then the same conclusion can be drawn for the entire sequence. □ 

Remark 1. In some cases value e~ 1 X e doesn't actually depend on e. The latter takes place if 
and only if the three-points problem 

(K(x)u')' + fiR(x)u = for x G I a U I b , u{a) = u(b) = u'(-0) = u'(+0) = (23) 
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has an eigenfunction u that is continuous at x = (for a certain eigenvalue /i). This situation is 
possible, for instance, in the case a = —b when there exists even eigenfunction of the Dirichlet 
problem on (—6,6). Then a trivial verification shows that X s = e\x is an eigenvalue of (J2J)- (JHI) 
with the eigenfunction u £ = u for all e G (0, 1]. 

Corollary 1. Restrictions of eigenfunction u £ j to the intervals I a and It, converge towards the 
corresponding restrictions of eigenfunction u in H 2 (a,0) and H 2 (0,b) respectively. 

Proof. Set u £ = u e j. We consider equation ([2]) in the form u" = —k'k~ l u' e — fi e rk^ 1 u e in I a . 
Then from Theorem [2] we have 

u" -> -k'k~ l u' - firk^u in L 2 (a,0), (24) 

where u is an eigenfunction of ( fl5l) . From (TT51) it follows that the limit ( 124"1) is exactly the 
second derivative of the limiting eigenfunction in I a . The proof for interval is the same. □ 



3. Formal Asymptotic Expansions of Eigenvalues and Eigenfunctions 

3.1. Asymptotics of Simple Eigenvalues. In this section we construct the complete as- 
ymptotic expansions of eigenvalues A £ and eigenfunctions u E . We begin with the examination 
of eigenvalues A £ for which the limit \x = lim \ e Je is a simple eigenvalue of operator A. Clearly, 

J e—>0 J 

H depends on j, which we do not indicate for the sake of notation simplicity. The asymptotic 
expansions of the eigenvalues and the corresponding eigenfunctions are represented by 

A £ ~e(/i + + ■■■ + e n v n + •••), (25) 

/ \ J Vo{ x ) + £ Vi( x ) H 1" £ n Vn(x) H for x G I a , , * 

I z {x) +ezi{x) ^ h e z n {x) H tor x G J&, 

where /x is an arbitrary eigenvalue of limit problem (IToT) . Then 

for x G / a , /oT\ 

=S fir ^ r ( 27 ) 




for x E h 

is the corresponding eigenfunction of ([151) as it follows from Th. [2j Since in this section we 
treat only the simple eigenvalues /i, according to Th. [1] we only consider here two possible 
situations: fi G 0"(Ai)\cr(A 2 ) and fi G a(A2)\cr(Ai) . 

3.1.1. Case /i G cr(yli)\o"(A2). We fix the corresponding eigenfunction ?/ of operator A\ such 
o 

that J ry 2 dx = 1 and yo(0) > 0- Since \i doesn't belong to the spectrum of A 2 there exists a 

a 

unique solution z Q to the problem 

{xz' o y + (ipz = in I b , zo(0) = y (Q), z Q (b) = 0. (28) 
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An easy computation shows that the next terms of the expansions are unique solutions to the 
recurrent sequence of problems 

n-l 

(W + (*ry n = -v n ryo -rJ2"j Vn-j in J , 

3=1 ( 29 J 

.y n (a) = 0, (ky' n )(0) = (xzUXO), f\y n y dx = 0, 



n 

'xz'J + ppz n = -p E *n-j in 

j=i (30) 

.*n(0)=J/„(0), z n (&) = 

with v n = — (x^_ 1 )(0)yo(0) for n = 1,2, ... . The last formula for u n is obtained as the 
solvability condition of (1291) . Note that all solutions y n , z n are smooth functions. 

Remark 2. It might happened that z' (0) = (cf. the proof of Th. [2]). In this case function u 
defined by (l2"Tj) is exactly an eigenfunction of the perturbed problem for each e G (0, 1]. Then 
the construction of asymptotics is interrupted and we can state that there exists an eigenvalue 
A e = ep for all e > 0. The corresponding eigenfunction 



u £ (x) 



y Q (x) for x G I a , 
Zq{x) for x G h 

doesn't depend on e. 

3.1.2. Case p G a(A 2 )\cr(Ai) . This situation immediately implies yo = (cf. the proof of 

b 

Th. [H part (ii)). We fix the corresponding eigenfunction z Q of A 2 such that J pz^dx = 1 and 

o 

z' (0) > 0. A trivial verification shows that the next terms of expansions ( 1261) are the unique 
smooth solutions to the problems 

n-l 

{ky'J + pry n = -rJ2 v i Vn~j in I a , 

3=1 

.y n (a) = 0, (V„)(0) = (x<^ 1 )(0), 



n-l 

>4)' + /Wn = -^nP^O ~ PH V 3 Z n-j hi /{ 

J'=l 



,^n(0) = y n (0), z n (b) = 0, j*pz n z dx = 0, 
with i> n = —(xz' )(0)y n (0) for n = 1, 2, . . . . Such choice of v n assures the solvability of ([31 



(31) 



3.2. Asymptotics of Double Eigenvalues. In this subsection we treat the case when for 
two successive eigenvalues A^ and A^ +1 the corresponding ratios £ -1 A^ and e~ 1 X e j +1 converge to 

the same limit p. It is obvious that p must belong to the intersection cr(Ai) U cr(A 2 ). Let us 
assume that the eigenvalues and the corresponding eigenf unctions admit expansions 

A £ ~ e(p + y/eux + ev 2 + ■ ••), (32) 



^ . y/ewi(x) + ew 2 (x) H for are (a, 0), 

v (x) + y/evi(x) + ev 2 (x) H for x G (0, 6), 
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because the eigenvectors of operator A that correspond to double eigenvalues \i have the form 
Vfj, = (0, i>o) (see Th. [[]). Substituting (132]) . (133]) into the perturbed problem we obtain 

(xv'q)' + fipv = in I b , v (0) = v (b) = 0, (34) 

{kw[y + firwi = in I a , wi(a) = = 0. (35) 

We fix /x G o~(A\) U o~(A 2 ) and introduce the functions 

U(x) = 1 ° for^e^a [/ ( x ) = j™*^ fo r ^e/ a 6 
1 f(x) for x E h ' )v*(x) for x £ h 

that correspond to the eigenvector and adjoined vector of ^4 (cf. vectors and in Th. [T]). 

Here w is an eigenfunction of ( 134]) such that J*„ 6 pt> 2 = 1, ?/(0) > and adjoined vector U* is 
chosen such that (U, UA L 2 (R,(a,b)) — 0. We also introduce an eigenfunction w of (1331) such that 
rw 2 dx = 1 and w;(0) > 0. It follows that i>o = aw an d = P w with certain constants a 
and /3. In addition, a must be different from zero. The next problems to solve are 

(*cv'i)' + npv\ = —v\apv in lb, fi(0) = /3w(0), = 0, (37) 

(kw' 2 )' + firw 2 = -ui/3rw in I a , w 2 (a) = 0, fc(0>4(0) = ax(0)i/(0). (38) 

In general case both problems (137j) and (1381 are unsolvable, since /x belongs to the spectra 
<r(y4i) and (7(^2) at one time. Hence we have to apply Fredholm's alternative for both the 
problems. After multiplying equations (138]) and (1371) by eigenfunctions v and w respectively 
and integrating by parts, one yields the common solvability condition: 

M R ' ( 39 ) 



u J \ p J ~ L \ (3 

where uj = (xwv')(0) is positive. Since the first component of vector 7 = (a.,0) must be 
different from zero, —V\ is an eigenvalue of the matrix in (1391) . Therefore if either V\ — ui and 
7 = (1,-1) or V\ = —uj and 7 = (1,1), then problems (1371) . (1381 admit solutions. Moreover, 
functions vivo* and V\v* solve problems (135]) and (1371) respectively for both values of v\. Actually 
these problems imply immediately (A — fi)U* = uU. In other words, the first corrector is an 
adjoined vector of A that corresponds to the eigenvector ujU. It causes no confusion that we 
use the same letters U, U* to designate a function of L 2 (a, b) and a vector in C 

Summarizing, we formally demonstrate that there exists a pair of closely adjacent eigenvalues 

and A^ +1 that admit the asymptotic expansions 

Xj = efi- e 3/2 uj + 0(e 2 ), = e/i + e 3/2 uj + 0{e 2 ), as e -> 0. 

As of asymptotics of eigenfunctions we have 

u £ j(x) = U(x) - s/euU*(x) + 0(e), u eJ+ i(x) = U(x) + y/eujU*(x) + 0(e). 

These eigenfunctions subtend an infinitely small angle in L 2 -space as e — * 0. Hence u e j and 
u £ j + i stick together at the limit. The latter gives rise to the loss of completeness of the limit 
eigenfunction system. 

Suppose that v\ = uj and 7 = (1,-1). Then we will denote by V\ and W 2 such solutions of 
the problems that J pV±vdx = and j\]V 2 w dx = 0. We see at once that —V\ and —W 2 are 
solutions of (137]) . (1381 for v\ = —uj and 7 = (1, 1). 
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From now on we distinct two branches of expansions (132]) 



A^ ~ e{n + ^uj + ev+ + ■ ■ ■ + e n ^ 2 v+ + ...), 



(40) 



and the corresponding branches of ( }33j) are 

±, , _ f TVew(x)± ew^{x) + --- + e n l 2 wt{x)..., x e I a , 

e y ' \ v(x) ± y/evf{x) + evf(x) H h £ n/2 i>*(:r) . . . , x e I b . y ' 

All coefficients are endowed with indexes + or — if they depend on the choice of the sign of the 
first corrector V\ = ±u. Note that the high order correctors in (]40p . (I4ip have to be calculated 
separately for both the branches. We now turn to the case v\ — u and find coefficients v£, 
and u+. To shorten notation, we omit upper index "+" for a while. Next, we see that 
problems (j57j) and (J38]) admit many solutions Ui = Vi + a.\V and W7 2 = W 2 + Piw, where ax, fix 
are constants. These constants can be obtained from the consistency of problems 

{xv' 2 )' + fxpv 2 = -vxp (Vx + axv) - v 2 pv, x E I b 
v 2 (0) = W 2 (0)+(3xw(0), v 2 (b)=0, 

(kw' 3 )' + /irw 3 = -v x r (W 2 + Pxw) - v 2 rw x , x e I a , , 

w 3 (a) = 0, k(0)w 3 (0) = x(0)(Vx + axv)'(0). 

The solvability conditions for problems (142]) and (143]) . which arrive from Fredholm's alternatives, 
can be represented as a linear algebraic system 



vx u \ ( a x \ _ ( (xW 2 v')(0) + v 2 
w vx \ Pi J V (*™Vx)'(0) - v 2 



(44) 



The system has solution if and only if v 2 = | (xwVx — xW 2 v') (0). After the solvability condi- 
tion is satisfied, system (144"]) has a partial solution ax = 0x — j- {xwV{ + >cW 2 v') (0) and prob- 

b o 

lems (1421) and (|43[) admit solutions V 2 and such that j pV 2 v dx — and J rM^to <ix = 0. 

a 

Therefore, all other solutions of (|42|) and (l4"3]) allow the representation v 2 = V 2 + «2 W and 
u; 3 = W 3 + p 2 w with real constants a 2 , (3 2 . 

We construct the general terms of expansions ( 1401) and ( 1411) as solutions to the problems 

n 

j=i ' (45) 
u n (0) = w n (0), = 0, 

n 

i=i (46) 
w n+1 (a) = 0, (&w n+1 )'(0) = (xf n _i)'(0), 

with 

^n-i = Ki-i + and u> re = W n + (3 n -xW, (47) 
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where V n -\ and W n are solutions of the previous problems chosen accordingly to the orthog- 

b o 

onality conditions J pV n ^\V dx = and J rW n wdx = 0, n > 2. Constants a n _i and j3 n _i we 

o a 
find from the solvability conditions for (1451) and (1461) given by 



/ n_1 \ 

' (KW n v') (0) + £ VjO^i + I/, 



n-1 



V 



(4* 



(xwK-l) (°) + E ~ V n 

j=2 / 

The latter has a solution if and only if v n — | (xw 

K-i - »W n i/) (0). Then system @E5 has 

n-1 

a partial solution a n _i = /5 n _i = ^ \KvSV' n _ x + xW n v'j (0) + - ^ u j a n-j- Substituting the 

constants into (I47p we finish the general step of recurrent algorithm. Hence, after coming back 
our natation we obtain all coefficients v£ and of series (14"U|) and (14"T1) . 

Similarly, we can construct the coefficients v~ , t>~ and u>~ of series (j4"U|) and (HIT) . Then, 
by induction we get that for any natural n the coefficients satisfy relations v~ = ( — l) n i/+, 

= (~l) n < and ^n = (-l) n ^n- 

4. Justification of Asymptotic Expansions 

Let C £ be he weighted L 2 -space with the scalar product and norm given by ([6]). We also 
introduce space TC £ as the Sobolev space Hq(cl, b) with scalar product and norm 



(M) e = k<f>'i>'dx + e xfi^'dx, U\\n e = VWWe- (49) 

J a JO 

It is easily seen that 

cU\\ < U\\e < Ce~^U\\, ce^Uh < U\\ ne < CU\\ U (50) 

where || • || and || • ||i are standard norms in £,2(0,, b) and Hq(ci, b) respectively. 

For the sake of completeness, we introduce here below the classical result on quasimodes. 
Let A be a self-adjoint operator in Hilbert space H with domain T>(A) and a > 0. 

Definition . We will say that pair (//, u) G K x T>(A) is a quasimode with accuracy to a for 
operator A if || (A — ^I)u\\h < a and \\u\\h = 1. 

Lemma 1 (Vishik and Lyusternik). Suppose that the spectrum of A is discrete. If (/j,,u) is a 
quasimode of A with accuracy to a, then interval \p — a, \i + a] contains an eigenvalue of A. 
Furthermore, if segment \ja — d, \x + d\, d > 0, contains one and only one eigenvalue A of A, 
then \\u — v\\h < 2d~ 1 a, where v is an eigenfunction of A with eigenvalue \, \\v\\h = 1. [HE] 

4.1. Simple Spectrum. We will denote by A e n = e (/x + ev\ + • • • + e n u n ) and 

y (x) +ey 1 (x) H h e n y n (x) for x G I a 

z (x) + ez\(x) + • • • + e n z n {x) for x G h 

the partial sums of series (1231) . (1261) . The perturbed problem is associated with self-adjoint 
operator A £ = -±£k £ £ in C £ with the domain V{A £ ) = {f E H: (kf')(—0) = e(x/')(+0)}, 
where coefficients k £ , r £ are given by ([I]) for m — 1. 



f4,n(^) 
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Theorem 3. If Hj G 0"(Al)\<7(A 2 ) ; then eigenfunction u £ j of ([S])-© with eigenvalue Xj con- 
verges in H 1 (a, b) towards the function 



u(x) 



y{x) for x G I a 
z(x) for x G h, 



where y is an eigenfunction of the problem {kg')' + /iry = in I a , y(a) = y'(0) = with 
eigenvalue \Xj, and z is a unique solution of the problem (xz')' + fij pz = in h, z(0) = y(0), 
z(b) = 0. 

If z'(0) = 0, then X £ j = efij and u £ j = u for all e > 0. Otherwise X £ - and u £ j admit 
asymptotics expansions (j25]L obtained in \3.1.1\ for \i — fij. Moreover, the estimates of 

remainder terms hold 

I^A* - (jjij + evi + • • • + e n v n )\ < c n e n+1 , (51) 
\\u £jj -$ £ U £>n \\ HHaib) <C n e n+ \ (52) 
where "d £ is a normalizing multiplier with strictly positive limit as e — > 0. 

Proof. The case z'(0) = was considered in Remarks [H and [2j Suppose that z'(0) ^ 0. We 
first check that the the series being constructed in 13.1.11 give us the quasimodes with accuracy 
to an arbitrary order. It follows from (1291) . (I30p that 

\r:\k e U' e J + K,nU £ ,n\ < C n E n + 2 (53) 

in [a, b] uniformly, U £>n (a) = U £j Jb) = 0, U £tTl (-0) = £/ e ,„(+0) and 

&.» = (kU' £ , n )(-0) - e(^,„)(+0) = 0(e n+1 ), e -> 0. (54) 

Note that U e>n doesn't belong to the domain of A £ since (3 £>n is different from zero in the 
general case. Set <f)(x) = x(- — 1) for x G (a, 0) and 0(x) = elsewhere. Then V e , n = U £>n + 
P £j n4> belongs to V(A £ ) and a simple computation gives || A:K,n — A e n V^ jn || e < c n e n+3 ^ 2 . Hence 
(A £jn , K )n /||V^ )n || e ) is a quasimode of operator A £ with accuracy to c n e n+2 because ||V e ,n||e = 
0{e~ 1 ^ 2 ). According to the Vishik-Lyusternik Lemma there exists an eigenvalue A e of A £ such 
that |A e — A e „| < c n e n+2 , which establishes floTl) . Moreover, there exists an unique eigenvalue 
A e = A^ with such asymptotics by Theorem [2J Next, for a certain d > segment [A e n — 
de,A e , n + de] contains one and only one eigenvalue of A £ . Repeated application of Lemma [1] 
enables us to write 1 1 1 1 1 1 ^ 1 ' u e ~ IlK^II^ 1 • K,n|| < 2c n d~ 1 e n+1 , where u £ = u £ j. Hence, by 



\\ u e\e w 
Ue ~ Try fP^e.n 



E.n £ 



<^\\u £ \\ £ s n+1 <C n s n+1 / 2 

d 



and $ £ = y converges to 1 by Theorem [2j 

Pair (A e , u £ ) satisfies identity (u £ , ip) £ = X e (u £ , ip) £ for all ip G Hq(o, b). Similarly, (V £jn , if)) e = 
Kn(Ve,n,ip)e + « e (^), where \a £ (ip)\ < ce n+1/2 \\ij\\ ne . The latter gives 

\\u £ - fi £ V £in \\ He < A £in ||w e - $ £ V £ , n \\ £ + |A £ - A e , n | ||lt e || e + \a £ (u £ - CK,n)| 

< 2 N C n e n+3 ' 2 + c n \\u £ \\ e n+3 ' 2 + ce n+l l 2 \\u £ - # £ V £>n \\ Hs 
and consequently \\u e — ^eV^nllwe < C n e n+3 ^ 2 . From this and (1501) we thus get estimate (1521 . □ 



ulx) 
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The same proof works for the rest part of the simple spectrum of A. 

Theorem 4. If fij G a (A 2 )\cr (Ai) , then eigenfunction u £ j of ©-([S]) with eigenvalue Xj con- 
verges towards function 

for x G I a , 
z(x) for x G h 

in H l (a,b), where z is an eigenfunction of the problem (xz 1 )' + fi pz = in 1^, z(0) = 0, 
z(a) =0 with eigenvalue fij. Moreover Xj and u £ j admit asymptotic expansions (1251) . ( 1261) 
obtained in \3.1.S\ for fi = fij with the estimates of remainder terms 

|£~ X A* - (nj + ev x + h e n v n )\ < c n e n+1 , \\u £d - $ £ U £ , n \\ H i {ayb) < C n e n+1 . 

Here d £ is a normalizing multiplier that converges to a positive constant as e — > 0. 

4.2. Double Spectrum. We introduce the partial sums of (140]) . (|4T!) 

\% n = e{nj ± e 1/2 u + ei/± + ■ • • + e n/2 v±), (55) 
TT ± / =R 1/2 ^ + H h £ n/2 w^ for x e I a , , 

U e,n | „ , i „l/2„,± , i _n/2„.± f_ „. ^ r ^"J 



£ ' n \v + e 1/2 vf + ■■■ + e n l 2 vt for x G J 6 
with all coefficients constructed in Section 13.21 for certain double eigenvalue \i = fij = fij+i- 

K.n and U t,n 



Set V^ n = Uf n + Pf n 4>, where (3 e n and (3+ n are residuals in condition (TJJ for U £ n and U £n 



respectively defined similarly as in fl54l . Moreover, (3f n = 0{e (jl+l ^ 2 ) as e — > 0. 
Analysis similar to that in the proof of Theorem [3] leads to the following result. 

Proposition 3. The pairs (A~ n , V^~ n /|| V^J| e ) and (A+ n , V^/||V£JJ| E ) are quasimodes of oper- 
ator A £ with accuracy to c n e a l 2 . 

Proposition 4. There exist two closely adjacent eigenvalues X~ and Xf of (EP-([3]) with the 
asymptotics 

^ = fij ± y/Iu + eu± + • • ■ + e n/2 v± + 0(e (n+1)/2 ), (57) 

where fij is a double eigenvalue of operator A and tu, u k were defined in Sec. \3.2l 

Proof. From Proposition [3] and the Vishik-Lyusternik Lemma it follows that there exists at 
least one eigenvalue of A £ in each e n l 2 - vicinity of Aj n and A+ n . Moreover, — A^J < c n e n l 2 . 
Evidently, eigenvalues A~, Xf are different, because A+„ — A~ n > ue 3 ^ 2 and e n l 2 - vicinities of 
Aj n and Af n don't intersect for n > 3 and sufficient small e. In fact, |A+ — A~| > ce 3 ! 2 for 
certain positive c. We conclude from \Xf — A^ n+3 | < c n+ ^ n+3 ^ 2 that 



A ± 



n+1 \ "V "+k . _r_ . n+1 



k=l 



which establishes f)57p . □ 

We consider two planes in L 2 (a,b). Let 7r be the root subspace that corresponds to double 
eigenvalue fii and n(e) be the linear span of two eigenf unctions u~ and uf that correspond to 
eigenvalues A~ and A+. These eigenfunctions as above are normalized by (l2lil . 
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Theorem 5. The root subspace ir is the limit position of plane ir(e) as e — > 0, namely \\P n (e) ~ 
P n \\ 0, where P n (e) an d P-k o,re the orthogonal projectors onto planes n(e) and n respectively. 

Proof. Nevertheless both eigenfunction u~ and u+ converge to the same limit being the eigen- 
function of A with eigenvalue p.,-, the n £ has regular asymptotic behaviour as e — > 0. We choose 
new L 2 (R, (a, 6))-orthogonal basis in n(e): f £ = \{u+ + u~), g £ = ^— ^ (uf — u~). 

By Theorem [2] the first vector f £ converges in L 2 towards eigenfunction U £ 7r given by (j36j) . 
Next, function g £ solves the problem 

/\/ Aj~ A_ — A_ - ■ -r / /\/ Aj~ A. — Aj~ - ■ T 

( k 9e)+ T rg e = ^— /fr u £ m/ a , + _ We = -^--^ m J 6> 

k (/ e (a) = 0, &(&) = (), ^(-0)=y e (+0), (A;^)(-0)=e(^)(+0). 

Since £ _1 A+ — > p.,-, e _3 / 2 (A+ — A~) — > 2lj by (137]) and the right-hand side is orthogonal to 
the eigenfunction u+ in £ e , one obtains that norms ||flf e ||ifa( 0) o) an d 11^11^(0,6) are bounded as 
£ — > 0. Taking into account Corollary [I] we can assert that each converging subsequence g £ > 
converges as e — > towards a solution of the problem 

Hj rg = in J a , (xg')' + p,,- p# = -pv in J 6) 
0, (7(6) = 0, g(-0)=g(+0) 7 ^(-0) = 0, 

because u~ converges to eigenfunction U, which equals v in 1^ and vanishes in I a . Hence, 
all partial limits of the second basis vector g e have to be the adjoined vectors corresponding 
to the eigenvalue Pj. In fact, by orthogonality of f £ and g £ these limits belong to the line 
{at/* | a £ R} C it , which is orthogonal to U (see ( 1361) for definition of [/*). □ 

Indeed, in previous statements A~ = A|, A+ = A^ +1 and u~ = u £ j, uf = u £ j + i, by Theorem 
[2j Next theorem summarizes all information on bifurcation of the double spectrum. 

Theorem 6. Let fij £ cr(Ai) ncr(A 2 ) be a double eigenvalue with eigenfunction U and adjoined 
function U* given by (|36|) . fij = fij+i- Then both eigenfunction u £t j and u e j + i converge to the 
same eigenfunction U and the difference -^(u £J+1 — u £ j) converges to adjoined function 7 {7* 

for certain 7 7^ 0. Besides, A7 = A|, A+ = A| +1 and u E j, u e , 3 -+i admit asymptotic expansions 
( HOI) . ( HIT) derived in Section (2H /or p = pj. 27ie estimates of remainder terms hold 

\8^\f - (p, ± v^a; + er^ + ■ ■ • + £ n/2 ^) | < c^ n+l ^\ (58) 

IKj - K U 7,n\\HHa,b) < C~£—, \\u Etj+1 - $t U tn\\m{a,b) < E~ , (59) 

where $f are normalizing multipliers with strictly positive limit as e — > 0. 

Proof. It remains to prove estimates (|59l) . From (l58l) and Theorem [2] it may be concluded that 
for certain d > and n > 2 interval [A~ n — cfe 2 , A~ n + de 2 } contains eigenvalue A^ only. In view 
of Prop. [3] and the Vishik-Lyusternik Lemma, we have 



II V £,n\\e 



< 2 -^\\u £ \\ £ e^<C n e^. 

n-4 



As in the proof of Theorem [3] we can obtain \\u e> j — , d~U £n \\H 1 (a,b) < C n e^~ . Since all the 
coefficients of sum U~ n are bounded in if 1 (a, 6), the first estimate (159]) follows from the last 
inequality with n replaced by n + 5. The same proof works for u £ j + i. □ 
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